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1 3.  ABSTRACT  (Modmum  200  wonls) 


Conventional  integrated  GPS/INS  navigation  systems  model  pseudo  range  (PR)  and  delta  range  (DR)  residuals  as 
functions  of  position,  velocity,  user  clock  bias  and  user  clock  drift  errors.  Any  dependence  of  PR  and  DR  residuals  to  atti¬ 
tude  errors,  due  to  differences  between  the  navigation  solution  reference  point  and  the  antenna  phase  center,  is  ignored. 
Attitude  errors  are  estimated  through  correlation  to,  primarily,  velocity  errors.  This  correlation  develops  due  to  the  cou¬ 
pling  of  the  attitude  errors  to  the  velocity  errors  via  the  presence  of  specific  forces. 

The  simplification  of  ignoring  any  dependence  of  residuals  to  attitude  errors  in  their  estimation  is  a  good  assumption 
for  most  terrestrial  applications  where  specific  forces  are  present  and  may  be  enhanced  through  appropriate  maneuvers. 
The  application  of  GPS  to  exoatmospheric,  free  falling  conditions,  however,  is  an  emerging  field.  Under  free  falling  condi¬ 
tions,  however,  the  specific  forces  are  identically  zero  and  thus  the  ability  to  estimate  attitude  errors  through  correlation  to 
velocity  errors  vanishes. 

This  paper  describes  the  development  of  an  approach  to  estimate  attitude  errors  in  vehicles  with  changing  attitude 
independent  of  the  value  of  the  specific  forces.  The  paper  includes  derivation  of  a  more  general  GPS  observation  matrix 
which  captures  any  attitude  and  attitude  rate  information  that  was  lost  in  the  conventional  derivations  and  simulation 
results.  The  technique  developed  here  has  the  advantage  of  being  readily  implementable  in  a  conventional  GPS  receiver 
with  only  minor  software  modifications. 

Published  in  Pmreedings  of  the  Forty-Eighth  Annual  Meeting  of  the  Institute  of  Navigation,  July  1992. 
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ABSTRACT 

Conventional  integrated  GPS/INS  navigation 
systems  model  pseudo  range  (PR)  and  delta  range  (DR) 
residuals  as  functions  of  position,  velocity,  user  clock  bias 
and  user  clock  drift  errors.  Any  dependence  of  PR  and 
DR  residuals  to  attitude  errors,  due  to  differences  between 
the  navigation  solution  reference  point  and  the  antenna 
phase  center,  is  ignored.  Attitude  errors  are  estimated 
through  correlation  to,  primarily,  velocity  errors.  This 
correlation  develops  due  to  the  coupling  of  the  attitude 
errors  to  the  velocity  errors  via  the  presence  of  specific 
forces. 

The  simplification  of  ignoring  any  dependence  of 
residuals  to  attitude  errors  in  their  estimation  is  a  good 
assumption  for  most  terrestrial  applications  where  specific 
forces  are  present  and  may  be  enhanced  through 
appropriate  maneuvers.  The  application  of  GPS  to 
exoatmospheric,  free  falling  conditions,  however,  is  an 
emerging  field.  Under  free  falling  conditions,  however, 
the  specific  forces  are  identically  zero  and  thus  the  ability 


to  estimate  attitude  errors  through  conntiation  to  velocity 
errors  vanishes. 

This  paper  describes  the  development  of  an 
apfBoach  to  estimate  attitude  errors  in  vehicles  with 
changing  attitude  independent  of  the  value  of  the  specific 
forces.  The  paper  includes  derivation  of  a  more  general 
GPS  observation  matrix  which  captures  any  ^titude  and 
attitude  rate  information  that  was  lost  in  the  conventional 
dcrivjuions  and  simulation  results.  The  technique 
develtqied  here  has  the  advantage  of  being  readily 
implemeniable  in  a  conventional  GPS  receiver  with  only 
minor  software  modifications. 

1.0  INTRODUCTION  AND  BACKGROUND 

In  integrated  GPS/INS  navigation  systems  the 
ebservation  matrix  is  conventionally  derived  under  the 
assumption  of  coincidence  of  the  antenna  phase  center  and 
the  reference  point  Consequently,  the  pseudo  range  (PR) 
and  delta  range  (DR)  measurement  residuals  are  modeled 
to  have  coupling  only  to  position,  velocity,  user  clock  bias 
and  drift  error  states.  Potential  coupling  to  attitude  errors 
or  attitude  rate  errors  is  ignored.  Instead,  residuals  are 
simply  corrected  for  the  antenna  lever  arm.  The  needed 
attitude  information  for  antenna  lever  arm  correction  is 
generally  provided  by  estimation  of  attitude  erron  and 
propagation  of  the  corrected  attitude  "total  state"  (typically 
a  direction  cosine  matrix  or  quaternion)  to  the  time(s)  of 
measurement  Attitude  error  estimation  for  such  an 
implementation  depends  entirely  on  correlation  between 
the  attitude  error  and  the  velocity  error.  Specifically, 
contribution  of  the  attitude  error  to  the  velocity  error  is 
due  to  the  product  of  the  attitude  error  vector  by  the  skew- 
symmetric  matrix  of  specific  forces  within  the  dynamics 
F  matrix 

This  approach  has  been  successfully  applied  to 
many  applications  such  as  icrreslrial  navigation.  There 
are,  however,  applications  where  the  correlation  between 
attitude  error  and  velocity  error  vanishes.  For  example, 
when  an  exoatmospheric  host  vehicle  is  frce-faJling.  the 


specific  forces  are  zero  for  all  practical  purposes.  Then, 
although  GPS  may  succeed  in  estimating  position  and 
velocity  errors  fairly  well,  attitude  estimation  is  not 
possible.  Moreover,  the  correction  for  the  antenna  lever 
arm  becomes  flawed  as  attitude  errors  grow  due  to 
instrumentation  errors  such  as  gyro  bias  and  scale  factor 
errors. 

Typically,  GPS-based  attitude  error  estimation 
methods  rely  on  receiving  carriers  from  multiple  GPS 
satellites  and  processing  the  single  difference  of  their 
phases.  A  minimum  of  3  non-coUinear  antennas  is 
required  in  most  methods  as  well  as  considerable  receiver 
hardware  design  features  and/or  special  software 
processing  beyond  the  typical  GPS/INS  designs. 

An  alternative  approach  to  altitude  estimation  is 
presented  in  this  paper.  We  seek  to  exploit  the  attitude 
information  inherently  present  in  the  PR  and  the  DR 
measurements  obtained  with  a  single  offset  GPS  antenna 
mount«}  on  a  platform  undergoing  attitude  changes.  The 
attitude  information  is  primarily  recovered  through 
processing  DR  measurements.  The  PR  measurement  is 
less  useful  because  of  the  generally  high  code  loop  noise. 
The  implementation  of  this  approach  does  not  require  any 
special  receiver  hardware  and  only  a  few  extra  terms  in 
the  standard  GPS/INS  Kalman  filter  observation  matrix 
H.  For  "good"  performance,  an  integrated  GPS/INS 
system  is  required.  The  INS  provides  for  propagation  of 
the  attitude  estimates  in  between  GPS  updates.  In  some 
specific  applications,  however,  GPS  alone  may  suffice. 

An  overview  of  the  contents  of  this  paper  is  as 
follows.  Sections  2,  3,  and  4  contain  independent 
mathematical  developments  of  this  approach  to  attitude 
estimation.  Section  2  provides  a  general  mathematical 
develc?)ment  of  the  observation  matrix  H  which  is 
applicable  to  any  changing  host  vehicle  attitude.  The  H 
matrix  is  developed  without  reference  to  any  of  the  other 
models  that  are  part  of  a  GPS/INS  Kalman  Filter.  Section 
3  provides  a  mathematical  development  of  the  observation 
matrix  H  which  is  also  generally  applicable  to  any 
changing  host  vehicle  attitude.  This  development, 
however,  is  based  on  a  state  space  formulation  and  as  a 
consequence  it  employs  and  thereby  relies  on  prior 
developments  -specifically,  the  state  transition  matrix  and 
the  plant  noise  model.  Section  3  also  provides  a  quick 
comparison  of  this  method  to  the  more  familiar 
"interferomeuic"  GPS-based  attitude  determination 
methods.  Section  4  provides  a  simpler  analysis  of  attitude 
estimation  but  limited  to  a  spinning  body  scenario.  As  a 
result  it  facilitates  gaining  certain  insights  to  and 
identifying  sensitivities  in  spinning  body  applications. 
Section  5  provides  simulation  results  for  a  spinning  body 
application  using,  however,  the  more  general  formulation 
of  Section  3.  Section  6  provides  a  summary. 


2.0  GENERAL  DERIVATION  OF  J  HE 

observation  matrix  H 

This  section  provides  a  general  rnatbemaucal 
development  of  the  observation  matrix  H  which  is 
applicable  to  any  changing  host  vehicle  aiutudc  We  first 
derive  H  matrix  row  vectors  corresponding  to  PR  residuals 
for  a  single  "offset"  GPS  antenna  and  then  proceeds  to 
derive  H  matrix  row  vectors  for  DR  residuals,  it  ls 
recognized  that  for  Jimilations  in  the  antenna  lever  arm 
length,  the  code  loop  noise  will  generally  overwhelm  the 
relation  of  PR  residuals  to  attitude  error.  The  DR 
residuals,  however,  being  based  on  the  very  low  norse 
carrier  tracking  may  well  afford  estimation  of  attitude  and 
attitude  rate  errors. 

It  is  assumed  that  the  LNS  can  measure  the 
antenna  rotation  with  sufficient  accuracy  so  that  the 
unwanted  phase  shifts  due  to  the  antenna  rotation  as  well 
as  other  antenna  phase  imperfections  can  be  compensated. 
The  problems  and  issues  associated  with  the  corresponding 
compensation  of  DR  residuals  and  antenna  calibrations 
were  not  investigated  here. 

Pseudo  Range  Measurement 

The  geometry  of  the  problem  is  shown  in  Figure 
1.  With  respect  to  some  convenient  coordinate  frame 
(e.g..  Earth  Centered  Earth  Fued  -  ECEF)  designated  as 
L,,  the  PR  between  the  antenna  phase  center  and  the  jth 
satellite  is  given  by 


Figure  1.  Problem  Geometry 

where  Vp^j  represents  measurement  noise.  As  bias  in  the 
satellite  clock  is  not  modeled,  it  is  not  included  in  (1), 
We  will  suppress  any  designation  for  our  selected 
convenient  reference  frame  Z,  but  maintain  su[>crscnpts 
for  other  specially  designated  frames,  r*  is  the  random 
process  vector  providing  the  location  of  the  GPS  antenna 
but  as  represented  in  L,.  r^  will  be  rapidly  changing  due 
to  the  changing  attitude  of  the  host  vehicle  However, 
when  represented  in  the  host  vehicle  Ixxly  frame  (as 
),  it  will  be  constant,  well  known,  and  relalcd  to 
through  the  diieciion  cosine  matrix  C®  which  is  nuxlcled 
as  a  random  pnxicss 

r  (t)  =  C  (tl  r; 


(2) 


In  the  approach  of  this  section  it  is  convenient  to 
disregard  PR  and  DR  measurement  noise.  The  spatial 
differential  which  we  will  lake  with  respect  to  the  nominal 
or  predicted  value  (i.e.  evaluated  along  the  nominal 
trajectory)  corresponds  to  the  PR  measurement  residual. 
That  is, 


52j  =  5[eJ(r^*r^-r^)  1  ♦ 


(3) 


5b„ 


As  e,  is  of  constant  (unit)  magnitude,  Scj  for  small 
magnitudes  must  be  orthogonal  to  e^  and,  therefore,  to  r, 
+  Ia  *  I»i  to  which  e,  is  parallel.  Hence, 
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as  is  a  constant  (column)  vector.  From  the 

development  in  [1]  we  have 
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We  will  maintain  a  first  order  approxirnauun  to  the  spattal 
and  temporal  difference  pans  in  order  to  use  a  linear 
Kalman  filter.  Taking  the  tune  dtffereniiai  of  die  PR 
measurement  residual  we  have 

d{6z,)  =  die’  5r  )  -  d{e'  R  Sbi  ^  dcSB^i  <*) 
and  on  expanding 

dsSijJ  •  (9) 

-(daM  -  diSS/i 

By  lineanty  we  can  exchange  d  { •  )  and 
5  ( ■  )  for  all  except  6^  which  is  not  separable  with 
any  proper  meaning.  Since  ihc  rate  of  atuiude  error  is 
equal  to  the  angular  velocity  error 

D  ( 5d )  ■»  0)  ,  -  0)  =  5uj  ^  (10) 

d  (  56)  =5(5)  dt 
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where  R  is  the  skew-syimnetric  matrix  form  of 
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The  H  matrix  for  PR  measurements  only  but  referred  to 
the  GPS  antenna  phase  center  is  given  by 
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If,  instead  of  bringing  the  H  matrix  to  the  antenna 
phase  center,  had  we  derived  an  H  matrix  based  on  the 
antenna  phase  center  coinciding  with  the  reference  point, 
that  is  R  =  0  (but  still  using  measurement  residuals 
based  on  ^tenna  lever  arm  corrections),  the  resulting  H 
matrix  would  show  no  coupling  between  the  residuals 
and  6^  .  By  eliminating  then  the  column  in  (7) 
we  readily  obtain  the  usu5  H  matrix  encountered  when 
PR  measurements  are  presented  to  a  Kalman  filler 
implemented  in  ECEF.  The  more  general  H  mairi-x  given 
in  (7)  appears  to  have  promise  for  attitude  observability 
but  would  require  a  large  offset  distance  of  the  antenna 
phase  center  from  the  reference  point  of  estimated 
quantities  to  overcome  the  noise  in  the  GPS  receiver  code 
loops 


Delta  Range  Measurements 

For  a  DR  measurements  model,  we  seek  tlie 
spatial  difference  process  of  a  temporal  difference  process. 
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where  we  have  used 
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Integrating  from  t,-T  to  t,  under  the  assumption  of 
temporally  "well  behaved"  5r,,  5^  6v,,  and  6f,  and 
using  the  Theorem  of  Mean  Value  for  integration  we 
obtain  (sec  [1]  for  details) 
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where  C;-  T  £  t,,  tj,,  C^.,  t  ,  t,  <  c. 

The  weakness  of  this  method  lies  in  the  inability  to 
precisely  determine  the  appropriate  values  for  these  limes 
in  the  interval  (t,,  t,-T).  Incorporating  both  PR  and  DR 
residuals,  the  resulting  H  mauix  is 
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If,  instead  of  bringing  the  H  matrix  to  the  antenna 
phase  center,  had  we  derived  an  H  matrix  based  O'  the 
antenna  phase  center  coinciding  with  the  reference  pomi. 
that  is  R  =0  {but  still  using  measurement  residuals 
based  on  antenna  lever  arm  corrcctiors),  the  resulting  H 
matrix  would  show  no  coupling  between  the  residuals 
84>  and  .  By  ellaiinaiing  the  50  column 

in  (25)  we  obtain  the  usual  H  matrix  encountered  when 
both  PR  and  DR  measurements  are  presented  to  a  Kalman 
filter  implemented  in  ECEF. 

When  Z,  is  a  local  level  frame  we  generally 
want  the  terrestrial  velocity  of  the  user  rather  than  the 
velocity  as  perceived  by  an  observer  fixed  to  Ig.  Using 
terrestrial  velocity  v,  the  observation  matrix  H  becomes 

53,  Sr. 

1  a 

l  3 
{  i 
1  0 
1  T 
0  T 

t  i 

0  T 

(15> 

Wc  note  that  5w  includes  the  effects  of  gyro 
bias  5o)„  ,  gyro  sc^e  factor  error,  etc. 

3.0  STATE  SPACE  FORMULATION 

In  developing  a  measurement  model,  the  method 
presented  in  section  2  has  taken  a  direct  approach.  The 
strength  of  such  an  approach  is  that  it  has  remained 
independent  of  any  other  analysis  and  thereby  has  been 
exploratory  in  seeking  coupling  between  PR  and  DR 
measurement  residuals  and  the  error  states.  The  results 
for  the  PR  residuals  need  no  further  elaboration  and  are 
useful  as  derived.  With  regard  to  DR  residuals,  however, 
a  weakness  is  that  it  has  not  been  established  where 
certain  terms  in  (13)  need  to  be  evaluated.  For  example, 
when  the  time  interval  T  over  which  the  temporal 
difference  is  taken  is  small  such  that  the  change  in  R 
is  small,  the  problem  is  minimal  and  any  reasonable 
choice  of  evaluation  time  will  provide  useful  results.  For 
greater  efficiency  in  coupling  to  the  attitude  errors, 
however,  we  wish  to  maximize  AR  .  The  choice  of 

where  to  best  evaluate  R  in  the  mterval  (t,-T,t,)  is 
uncertain. 

In  this  section  we  utilize  a  presumed  existing 
dynamic  model  to  develop  the  temporal  differences 
required  for  the  relations  through  backward  transition 
matrices,  'fhe  strength  of  this  approach  is  that  there  is  no 
uncertainty  as  to  where  to  evaluate  the  terms  m  tlie 
resulting  relations.  The  weakness,  however,  is  dial  it 
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depends  on  a  previously  developed  dynamics  model  and 
is  subject  to  any  inadequacies  ui  this  model.  In  parucuku. 
uncertainly  in  the  backward  time  propagation  coninbutcs 
to  uncertainty  in  the  measurement  model. 

The  touil  state  dynamic  model  is  taken  to  be 
continuous  and  given  by 

x!tl  ♦Gttjwtt)  (16) 

where  w(i)  is  a  zero-mean  white  Gaussian  noise  process 
For  the  purpose  of  cither  a  linearized  or  extended  Kalman 
filter,  we  assume  a  nominal  trajectory  of  similar  form 


X  (c)  »  E(x  ICl.t]  (17) 

— "W  “"“tt 

The  spatial  differenc*  process 

5x(t)  =  x(t)  -  x,ltl  (!«) 

has  the  first  order  approximation 

5x(t)  =  Ftt;xjc)  I5x(t)  *  G(t)w(t)  (19) 
where 


,  ,  (x, c) 

F(t;x  (C)  ]  - 

•“  "  dx 


The  discrete  time  measurements  are  taken  to  be  a 
nonlinear  function  of  the  state  with  an  additive  white  noise 
sequent*. 


z.(ti)  ♦  v(ti)  (20) 

Associated  with  the  nominal  trajectory  is  the  sequence  of 
nominal  measurements 


z  (CJ  h[x  (tj  ,tj  (21) 


The  spatial  difference  process 


5i(6.)  ”  2.(6.) -.L,!-,)  *  .tl  -  h(x^  (C.)  .  t ;  1 -vC.,; 

(22) 

has  the  first  order  approximation 


5,z(t,)  =  H[t..;><^(t,)  ]5x(t,)  V(t.) 

where  t  x  ( c.  )  ,  ^  i 


dx 


(23) 

(24) 


As  a  discrete  time  proce.ss,  the  measured  DR 
corrupted  by  measurement  noise  Ls 

Zo,,(tJ  =R.(r,)  -  R,(t,-T)  *  (25) 

where 

I  *  3Jt,) 

The  nominal  DR  measurement  corresponding  to  the 
nominal  trajectory  is 

^  R.JtT  -  R,Jt.,-T)  (26) 

The  DR  mcasuremciK  residual  discrete  time  prcKcss  is 
then 


•  R,(C,)  -  R,.IC,)  -lR,(e,-T)  -R,.U,-Tt|  •  v... 

(27) 

or  from  Figure  I 

6i„,(C,)  -ir^(E,)  -Sr^lC,)  I  'SS.U,)  (28) 

il£,  I  *l£,  -5r^it,-T)  -5r^!c,-Tl  1  *5b.II,-TI) 

Using  the  skew-symmetric  form  R  { t )  for 
vector  r  { t )  ,  we  have 

. "An 

6r  (t)  =  -  R  (t)56(C)  (29) 

We  then  have  for  the  discrete  time  measurement  residual 
process 

-|r,(t,)-5r^tt,)*  R^(tj5*(t.l|  -  Sa.lc.) 


-(eMt.jR  (t.) -eMt  -  Ti  R  (t,-T;lK 

— I**  ‘  -nH*  *  — 

*£j.  7-)£Fit,,x^  it,}  )T 

-e,^,  ( t  i  -  T)  it,-  T)  K.  F  1 1 , ;  i  £  , )  )  T 
!T 

The  DR  measurement  results  of  this  method  were 
compared  to  those  of  the  first  method  in  detail  in  [!}.  The 
results  were  shown  to  be  identical  other  than  the  time  of 
evaluation  is  quite  clear  for  the  second  method.  This 
overcomes  the  major  weakness  in  the  first  approach.  The 
second  method,  as  previously  noted,  inaeases  uncertainty 
in  the  measurement  model  for  any  inadequacies  in  the 
dynamics  model  used.  Both  methods  may  be  extended 
beyond  a  first  order  model. 

Comparison  To  Standard  Attitude  Measurement 
Techniques 


-{l£,ttrT)l  -lz,(t.-TI-5l,(t,-T)  ♦R_^(t,-T)5l(C.-Tll  -Sajc,-T)} 

(30) 

We  define  3  x  n  matrices  K  and  K  and  row  vector 

sap 

(1  X  n  matrix)  such  that  we  can  relate  to  the  full 
error  state  vector  at  measurement  lime  t,. 

5r^(tj)  =  l^5x(ti) 

5i(Ci)  =F8x(t,) 

8B„(ti)  =  kjSx(tj> 

At  measurement  time  ti-T  we  obtain  the  corresponding 
quantities  but  using  the  backward  transition  matrix 

8r^(ti-T)  =  i^[<t'(ti-T,Ci)8x(ti)  -► 

S*(ti-T)  =  j^[<[>(Ci-T,  tj)5x{Cj)  ■••w„(Ci-T)] 

5b^(Cj-T)  =  k^Olti-T.tjlSxtti)  -w^(Cj-T)] 

(32) 

where  WoCij-T)  is  the  driven  response  at  vT  due  to  the 
presence  of  the  white  noise  in  (16)  during  the  backward 
interval  from  ij  to  tj-T.  ^(tf-T),  by  virtue  of  the  white 
noise  in  the  continuous  model,  is  a  white  noise  sequence. 
It  is  here  that  uncertainty  in  the  dynamics  model  is 
introduced  in  the  measurement  model. 

The  discrete  measurement  residual  process  is  now 

[R^-R^lt,))^]«x(c,)l  (33) 

-|lr,(C,-Tl|  -|r,lt,-TS-  R^IC,-T)£^] 

-kj  (  •  s,ie.-Tl)  • 

If  we  approximate  the  backward  transition  matrix 
by  the  first  order  form 

<P(t,-T,t,)  =I  -  F[ti;x  (tj)  iT  (34) 

••  ‘  ‘  •"n.n  '  — "  ‘ 

we  obtain  the  jth  row  of  H  corresponding  to  the 
measurement  residual  for  the  jth  satellite  as 


Although  a  comparison  to  an  exhaustive  list  of 
other  GPS -based  attitude  determination  methods  is  beyond 
the  scope  of  this  report,  a  quick  comparison  of  the  method 
described  in  this  report  to  "interferometric"  attitude 
determination  methods  (sec  [2]  and  the  references  cited 
there)  may  be  useful  at  this  point.  We  can  construct  a 
configuration  which  affords  a  convenient  general 
comparison  between  the  technique  presented  here  and  the 
more  familiar  attitude  measurement  techniques.  Fcr 
purposes  of  comparison  we  take  any  two  non-coincident 
GPS  antennas  of  the  "standard"  interferometric  method 
and  the  set  of  their  phase  differences  for  each  of  four  GPS 
satellites  to  comprise  a  spalial  GPS  Interferometer.  A  set 
of  DR  measurements  taken  over  the  same  time  interval  for 
four  GPS  satellites  using  a  single  antenna  in  effect 
provides  a  temporal  version  of  a  GPS  Interferometer  as 
defined  above.  While  more  restrictive  than  need  be  due 
to  use  of  recunive  cstimatiOT  techniques  which  preserve 
information  by  propagating  it  forward,  the  typical 
implementation  of  a  conventional  five  channel  GPS 
receiver  having  PR  and  DR  measurements  would  in  fact 
conform  to  such  a  configuration. 

A  brief  comparison  between  the  technique 
described  in  this  report  (denoted  as  "DR-Based  method") 
and  standard  GPS  interferometric  methods  is  provided  in 
Table  1. 

4.0  SIMPLIFIED  DR  ANALYSIS  FOR  A 

SPINNING  BODY 

In  this  section  we  investigate  the  observability 
afforded  to  attitude  errors  by  the  DR  measurement 
specifically  for  a  spinning  body  scenario  Restricting  the 
analysis  to  a  simple  spinning  body  scenario  facilitates 
gaining  certain  insights  to  and  identifying  sensitiviues  in 
spinning  body  applications.  Ihc  analysis  is  limited  to  the 


Table  1 


Cotvparison  of  Standard  Interferometric  GPS-based  attitude 
measurement  techniques  to  DR-Based  method 


Standard  (Multiple  Antennas)  j 

DR-Based  (Single  Antenna) 

Requires  multiple  antennas  and 
a  special  GPS  receiver  design. 
No  requirement  on  host  vehicle 
motion  relative  to  GPS 
satellite  constellation. 

Requires  only  a  single 
antenna  and  a  standard  CPS 
receiver  having  PR  and  DR 
measurement  capability. 

Requires  changing  attitude  of 
host  vehicle  relative  to  GPS 
satellite  constellation. 
Integrated  GPS/INS  generally 
needed  for  good  performance. 

Any  two  non-coincident 
antennas  receiving  carriers  of 
four  GPS  satellites  and  the 
single  difference  of  their 
phases  constitute  a  spacial 

GPS  Interferometer. 

Each  set  of  DR  measurements 
for  four  GPS  satellites  using 
a  single  antenna  on  host 
vehicle  undergoing  attitude 
changes  is  equivalent  to 
constructing  temporally  a  GPS 
Interferometer  havi.ng  as 
antennas  the  start  and  stop 
antenna  locations.  The 
accumulated  phases  of  the  DR 
measurement  set  are 
equivalent  to  the  set  of 
differences  of  received 
phases  at  the  start  and  stop 
antenna  locations. 

A  minimum  of  two  non-parallel 
GPS  Interferometers  (a  minimum 
of  three  non-collinear 
antennas)  is  required  for 
attitude  determination  at  any 
time. 

A  minimum  of  two  non-parallel 
equivalent  GPS 
interferometers  (two  DR 
measurement  sets)  required 
for  attitude  determination 
with  INS  propagation  of 
information  over  the  DR  time 
period.  Attitude  is 
determined  for  end  of  period 
and  is  propagated  forward  by 
INS. 

The  satellite  phase  errors  are 
removed  by  the  single  spacial 
difference  in  phase  between 
the  two  antennas  of  the  GPS 
Interferometer . 

Removal  of  satellite  phase 
errors  is  inherent  to  each  DR 
measurement . 

Ambiguity  resolution  for 
integer  cycles  must  be 
accon^jlished  for  the  single 
phase  differences.  Use  of 
redundancy,  search  procedures, 
triple  differences,  etc.,  are 
required. 

No  ambiguities  in  each  DR 
measurement  provided  no  loss 
of  carrier  lock  or  cycle  slip 
in  GPS  receiver  channel 
associated  with  GPS 
satellite . 

May  require  calibration  and 
compensation  for  antenna 
induced  error  due  to  phase 
pattern  differences  between 
the  two  antennas  constituting 
the  GPS  interferometer. 

Fundamental  phase  change  due 
to  attitude  changes  in  single 
antenna  must  be  compensated 
for.  Antenna  must  also  be 
calibrated  for  phase  pattern 
anomalies . 

Double  difference  (difference 
between  two  independent  single 
differences)  required  to 
remove  receiver  timing  errors, 
electrical  path  bias  errors, 
etc . 

Use  of  single  antenna  and 
standard  DR  measurement 
obviate  the  need  for  double 
differences . 

Direct  solution  of  attitude 
performed . 

Recursive  estimation  (using 
Kalman  Filter)  of  attitude 
errors  relative  to  nominal 
trajectory'  established  by 

INS  and  its  corrections 
(resets)  is  performed. 

case  of  an  initial  (constant)  attitude  error  in  a  rotating  but 
not  translating  vebicle. 

From  (15),  the  attitude  error  ^  is  related  to  the 
DR  residual  for  the  jth  satellite  via 

H  .  =  -(Ae^  R  *  AR  )  (36) 

— js*  — — jn  >“\n 

For  our  choice  of  reference  frame  and  short 
periods  of  time  the  change  in  the  unit  vector  e^  is  small 
so,  we  can  assume  that  Ae^  is  approximately  zero.  Then 
(36)  reduces  to 

— -jn 

The  contribution  of  the  attitude  error  ^  to  the 
DR  measurement  residual  is  obtained  by  multiplying  both 
sides  of  (36)  by  5^. 

H...56  =  AR  56  (38) 


Since  we  have  assumed  diat  L  is  snuill  compared 
to  r^.  is  approximately  e^ijual  to  e,  Finally,  assuming  we 
are  only  interested  in  the  y-component  of  die  atuiude 
error,  we  obtain  from  (37) 

orr  T 

«  2L  sint-j-)  (ej;,cosu:tj-.^)  ♦  e^..sina){  e  i  ] 

(41) 

From  (41)  we  conclude  that  h^  is  zero  if  any  of 
the  following  three  conditions  occur 

L  =  0  (42) 

ilE  =  NJt  (N  =  0,1,2 _ }  (43) 


e,.,, cos  (cot, -2^)  ♦  e^j^sin{ci)t,--2-i  =  0  W) 


From  (38)  we  conclude  that  the  contribution  of 
the  attitude  error  86  to  the  DR  residual  is  maximum  when 
the  vectors  Arj^,  and  ^  are  mutually  orthogonal.  The 
contribution  is  zero  if  any  two  are  in  the  same  plane. 
From  Figure  2  we  conclude  that  the  skew-symmetric 
matrix  is 

0  -Lsinot^  0 

§  =  Lsinwt^  0  -Lcosot,  (39) 

An 

0  Lcos<ac^  0 

j*t»  Vi 


Obviously,  the  condition  in  (42)  can  occur  if  and 
only  if  the  lever  arm  is  zero.  The  condition  in  (43)  can 
occur  if  there  is  no  rotation  or  if  the  rotation  rate  is  such 
that  a  multiple  of  360  degrees  is  swept  during  the  DR 
integration  time  T.  The  condition  in  (44)  can  occur  if 

e.j.  =  0  (4S) 

or 

(46) 

The  condition  in  (45)  occurs  if  the  satellite  is 
directly  above  (or  below)  the  vehicle;  i.e.,  if  the  vector  r^ 
in  Figure  2  has  only  a  Y  component.  The  condition  in 
(46)  occurs  if  the  angle  swept  by  the  antenna  lever  arm 
during  the  DR  integration  interval  is  symmetrically 
distributed  about  the  line  bb’  in  the  X-Z  plane  -  see  Figure 
2. 

We  can  express  (41)  in  tenns  of  the  inner  product 
between  and  the  antenna  location  vector  at  the  mid 
point  of  the  DR  integration,  r^(t,-T/2). 


Figure  2.  Geometry  of  Simplified  Example 

and  thus,  after  some  manipulation  of  trigonometric  terms, 
we  have 

0  --ZLsin  —  COSU  (ti-.| )  0 

”  -syra  0  ZLsin—sinu  (Cj-2) 

0  -sym  0 

(40) 

where  "sym"  denotes  an  entry  equal  to  the  symmetric 
entry  in  the  matrix. 


=  2  sin(i^)  l^j  (47) 

=  2  L  sin(i^)cos(p[e^j,X,  1  ) 

where  r,^(t,-T/2))  is  the  angle  between  the  unit  vector 
^  and  the  antenna  location  vector  r,^(tj-T/2). 

From  (47)  we  conclude  that  h^  is  sinusoidal  and 
has  a  maximum  peak  value  of  2L.  Thiis  value  is  obtained 
when  the  vectors  and  r*(t,-T/2)  are  parallel  and  the 
antenna  rotates  180  degrees  starting  and  ending  on  line 
bb’,  hjs,,  is  zero  if  the  vectors  e,^,  and  r.^(t,'T/2)  arc 
orthogonal. 


Suppose  now  an  Information  Kalman  Filter  is 
used  to  eslunatc  the  constant  altitude  error  using  DR 
measurements.  The  variance  of  the  estimation  error  at  4 
is  (see  [3]) 

i»0  Obn 

where  Cor  denotes  the  sigma  value  of  the  observation 
noise.  For  sufficiently  large  t^,  we  can  assume  that  the 
sum  can  be  evaluated  using  an  average  (and  therefore 
constant)  value  in  place  of  h^j^fij).  Then  assuming  the 
average  value  to  be  equal  to  half  the  peak  value,  we  can 
show  that  for  tt»0,  (see  [1]  for  details) 

{to  ^  ^ -  m 

L  sin{i^)  cosp 

We  note  that  (49)  is  only  a  bound  derived  under 
very  specific  and  limiting  assumptions.  Consequently,  it 
is  not  very  useful  in  detennining  attitude  uncertainty 
bounds  in  more  general  situations.  It  does,  however,  help 
in  revealing  the  general  dependence  of  the  attitude 
uncertainty  to  L  and  coT.  The  following  conclusions  can 
now  be  drawn; 

1)  The  uncertainty  in  estimating  attitude  errors 
with  an  offset  antenna  is  proportional  to  the  DR 
measurement  noise,  i.e.,  Cdr.  (Odr  ‘S  inversely 
proportional  to  the  SNR  at  the  input  of  the  carrier 
tracking  loop). 

2)  is  inversely  proportional  to  the  length  of  the 
antenna  lever  arm  L. 

3)  is  minimized  if  the  angle  swept  during  the 
DR  integration  is  180  degrees. 

4)  is  minimized  if  the  satellite  is  in  the  plane 
defined  by  the  rotating  antenna  and  the  angle 
swept  during  the  DR  integration  is  symmetrically 
distributed  about  the  vector 

5.0  SIMULATION  RESULTS 

This  section  summarizes  covariance  simulation 
analysis  results  for  a  spinning  body  using  the  more  general 
results  of  section  3.  Detailed  discussion  of  this  analysis 
is  given  in  [1],  The  scenario  consists  of  a  vehicle  moving 
in  an  exoatmospheric  trajectory  (free  fall)  and  spinning  at 
a  constant  rate  to  about  its  Y  body  frame  axis.  Figure  3 
shows  the  vehicle's  altitude  profile.  The  vehicle  is 
equipped  with  an  integrated  GPS/Sirapdown  Inertial 
navigator.  The  GPS  antenna  is  mounted  on  the  perimeter 
of  the  vehicle  at  a  distance  L  from  the  axis.  Therefore, 
the  antenna  rotates  at  the  vehicle  spin  rate  to.  To  an 


m 


Fijure  }.  Altitude  vs.  Time 

approximation,  the  scenario  here  resembles  the  condiutms 
assumed  in  ibe  analysis  of  section  4.  Consequently,  the 
insights  gained  in  section  4  carry  over  to  the  more  realistic 
conditions  of  this  section.  Typical  errors  were  assumed 
for  both  the  GPS  and  the  INS  (see  [I]  for  details). 

Results  were  obtained  for  a  spin  rate  of  co=450 
deg/sec  and  a  DR  integration  time  T=1  sec.  Assuming 
first  that  the  antenna  was  mounted  on  the  spin  axis  (i.e. 
L=0),  the  RSS  attitude  error  was  shown  to  grow  to  about 
30  degrees  as  a  result  of  a  50  ppm  scale  factor  enor  (see 
Figure  4).  The  inability  to  estimate  this  attitude  error  was 
the  direct  consequence  of  zero  specific  fences  in  free  fall 
and  therefore  decoupling  of  attitude  errors  from  velocity 
errors.  When,  under  the  same  conditions,  the  antenna 
lever  arm  was  set  to  L=7  inches,  the  RSS  attitude  error 
was  estimated  to  0.2  degrees.  Furibennore,  the  scale 
factor  error  was  now  estimated  to  within  a  fraction  of  a 
ppm  (see  Figure  5). 

Figures  6  and  7  show  the  sensitivity  of  the  final 
attitude  enor  to  the  lever  arm  length  and  spin  rate, 
respectively,  for  T=0.78  sec.  These  results  agree  with  the 
predictions  in  Section  4  (i.e.  an  inverse  dependence  on  L 
and  cosecant  dependence  on  to). 


Figure  4  RSS  Atuiude  L'rror  Profde 
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Figure  5.  Y-Axis  Gyro  Scale  Error  Profile 


Altitude  Error  (dog) 


Figure  6.  Final  RSS  Aionide  Error  Sensitivity  to 

Antenna  Lever  Arm  (T=0.78  sec,  L=7  in) 


The  method  presented  here  appears  promising  as 
a  GPS-based  attitude  estimation  technique.  It  imposes 
minor  implementation  impacts  but  requires  host  vehicle 
attitude  changes.  Consequently,  its  main  utility  is  in  space 
applications.  For  terrestrial  applications  some  benefit  may 
be  derived  when  rotary  motion  is  inherently  presenq  e  g., 
rotary  wing  aircraft 
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